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2.6 Oseen's improvement for slow flow past a cylinder 

2.6.1 Oseen's criticism of Stokes' approximation 

Is Stokes approximtion always good for small Reynolds number? 
The exact Navier Stokes equations read 

V-g = (2.6.1) 

q. Vg = -Vp + z/VV- (2.6.2) 
Let us estimate the error of Stokes approximation where the velocity field is given by 

= (l + 1^ - I) . <!• = -^Vs^0 (l - 1^ - I) . (2.6.3) 

In the far field r/a^ 1, the viscous stress is dominated by the last term 



The inertia term is dominated by 



W— ~ O 
dz \ 



(r 

— ^ - O - 



Hence the error is their ratio 

which becomes unbounded for rja ^ 1. Thus inertia cannot be ignored in the far field. 

The difficulty is severe for the two-dimensional fiow past a cylinder. By taking the curl, 
Stokes equation gives 

= 0. 

Since the body is a source of vorticity, (f would become unbounded logarithmically for large 
rja. This is certainly unphysical and is known as Stokes' paradox. 
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In view of the increasing importance of inertia in the far field, Oseen suggests that the 
hnear approximation of the inertia term, which is of dominant importance in the far field, 

— * — * 

and not in the near field, be added. Let q — Wk + q' the Oseen equations are: 

V • = (2.6.4) 

W^^-^ + iyV^q'. (2.6.5) 
oz p 

The added linear inertia term is of dominant importance in the far field and not in the near 
field. 

We shall demonstrate the Oseen approximation for the two-dimensional case of a circular 
cylinder. 

2.6.2 Oseen's theory for a circular cylinder 

Because of (2.6.4), the pressure is still harmonic 

V^p = 0. (2.6.6) 

Now the velocity can be expressed as the sum of a potential part V0 and a solenoidal part 
g", where 

q'^W(f) + ^' (2.6.7) 

It follows that 

VV = (2.6.8) 
and vorticity is only associated with q". Substituting Eq. () into Eq. () and setting 

P = Poo = -pW^ (2.6.9) 



we get 

w 

dz 

which implies 



dz 



= uV^^', (2.6.10) 



= uV'C (2.6.11) 
where ( = V x q". To solve for q" let us introduce a vector potential ka 

(^Vx{ka). (2.6.12) 

Now 

V X C = V X (V X A^a) = V(V • {ka)) - kV^a = V - WV. 

oz 



On the other hand 



V X C = V X (V X g^O = V(V • - V^g^' = -V^g^' = - 



W dq" 
u dz 



after using (2.6.10). It follows by equating the two preceding equations, 

W dq^' _ da -2 
v az az 

Now (2.6.12) does not define cr uniquely. Let us impose a condition 

.da 



Then we can integrate to get 
hence, 



az 



^ 7, V. ^ r 
q" = Va ka, 

V V 



q" — ka — 



In polar coordinates the velocity components are 

V da 



((l = a cos 9 



W dr 



II ■ a 9a 



The total disturbance is 



I 9(f) u da 

,^^_ + ,eos^--- 

, 1 90 ^ u 1 da 

It is convenient to rearrange (2.6.14) by introducing 

a(y, z) = e^^/^'^a(y, z) 

so that 



V'a 



We now introduce dimensionless variables 

(t) = Wa^, a = WT., a = Wt., q' = WQ, r = aR. 
Then the governing equatios are 

= 



(2.6.13) 

(2.6.14) 

(2.6.15) 
(2.6.16) 

(2.6.17) 
(2.6.18) 

(2.6.19) 
(2.6.20) 

(2.6.21) 
(2.6.22) 

(2.6.23) 
(2.6.24) 
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and 

(V^-e^)$ = (2.6.25) 
The velocity components are 

a$ , 1 a$ 

1 (9$ 119$ 

The boundary conditions are 

g^j + cos ^ = R = l (2.6.27) 

Qe- sin 9 = R=l. (2.6.28) 
The general solution to (2.6.24) and (2.6.25) is 

< , . cos 9 

^ = AolnR + Ai—— (2.6.29) 
R 

and 

a^CoKoieR) (2.6.30) 

We shall determine the coefficients AqAi and Cq approximately for small € = Re/2. 
On the cylinder it! = 1, we use the approximation for eR <^ 1 

2d2 



KoieR) - - (<7 + ^ In ^) /o(6i?) + + 



- + - In — ) +... (2.6.31) 



2 2 y 

(7 ^ Coe^-^"<^^^Xo(ei?) 

1 e/? 

^ -Co(l + ei? cos e + ---)(cr + - In — + ■■■) (2.6.32) 
Details of calculation for small eR are given for convenience 

9$ >il n 

dR^R-R-'''^ 

(aR \ 
(T + In — j cos 9 

-a sin 6' |i= +Co(l + ei? cos 6' H ) (^(t + In sin 6* 
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}_}_da_ 



+ 



2e 
1 



Co e cos 6^(1 + ei? cos 6*) ( 0- + In — ) + • 



+— Co(l + eR cos 61 + 



-5' + 



+ 



1 1 

2^R 



Co{-eR sin 6')(1 + e cos 9) (a + In 



On the cyolinder R — 1 condition on the radial velocity requires that 



(sR \ 
cr + In — j cos 9 

1 



-Co e cos 61(1 + e cos 61) ( cr + In - ) + 



-^Co(l + ei? cos 9) [-^) + 



Neglecting 0(e In e) and equating coefficients of constant terms an cos^ separately, 

^0 + 1^0 



-A, 



i-i)c„(....|) + f = -i 



Prom the condition Qg — sm 9, we get 



The final results are 



Co 



|-(<7 + ln f)- 



Ao = +- 



e|-(a + ln f) 



and 



Ai 



i-(<7 + ln|) 4 
As a physical deduction, we leave it as an exercise to show that 



drag force/length = y- 



(j — In 



so that the drag coefficient is 

Cd 



D 



2tx 1 
e i - a - In f 



(2.6.33) 
(2.6.34) 

(2.6.35) 

(2.6.36) 
(2.6.37) 

(2.6.38) 
(2.6.39) 
(2.6.40) 
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In the far wake ei? » 1 



^o(ei;)^,/^e-(l-J5 + ...) (2.6.41) 



the rotational part is 



Coe^^-^y^e-^ = ^oi/^e-^(^--^) (2.6.42) 



There is significant contribution only in the region ^ <^ 1. Since 1 — cos^ 

^ „-eR(l-cos6>) , Qo 



2 ' 



. C„^^e--/.,^.e„;^,-..«,^. (,,43) 

the wake is therefore parabolic in shape. In the far field, the potential part s emits fluid 
isotropically as a source at the discharge rate (27rAo). On the other hand the rotational wake 
is a fluid sink with the influx rate 



2Cn 



r xf^e--^y^Rde = ^ r e-^'d, = ^ = -2nAo (2.6.44) 
^0 V 2eR e Jo e 



in view of (2.6.33). Mass is conserved. 



